Propagators and phase structure of Nf = 2 and Nf = 2 + 1 QCD 
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We investigate the phase structure of QCD at finite temperature and chemical potential by solving 
a coupled set of truncated Dyson-Schwinger equations for the quark and gluon propagator. In 
contrast to previous calculations we take into account the full back-reaction of the quarks onto the 
Yang-Mills sector and we include the effects of strange quarks. We discuss the resulting thermal mass 
of the unquenched gluon propagator and extract order parameters for the chiral and deconfinement 
transition from the quarks. Our result for the temperature dependence of the quark condensate at 
zero chemical potential agrees well with corresponding lattice calculations. We determine the phase 
diagram at finite chemical potential and find a potential critical endpoint at (/j,f p ,T EP ) ~ (190, 100) 
MeV. 



I. INTRODUCTION 

The existence of a high temperature and/or den- 
sity phase with quarks and gluons as thermody- 
namically active degrees of freedom has been a ma- 
jor prediction since the early days of QCD. Over 
the years this topic has received a lot of attention 
from both, the theoretical and experimental sides. 
Today, high quality results are available from lat- 
tice simulations at zero chemical potential, see e.g. 
[l], [1[ and references therein. As it turned out, the 
low and high temperatures phases of QCD are not 
separated by a phase transition but rather contin- 
uously connected. By now, lattice QCD has firmly 
established the notion of a crossover at physical 
quark masses and zero chemical potential. 

Unfortunately the situation is much less clear at 
(real) finite chemical potential, where lattice calcu- 
lations are hampered by the notorious fcrmion sign 
problem. Although various extrapolation meth- 
ods agree with each other at rather small chemi- 
cal potential 0-0] > regions in the (T, /^)-plane with 
Hq/T > 1 are hardly accessible. Therefore, despite 
these efforts the basic structure of the phase dia- 
gram of QCD is not at all settled yet, see e.g. d,|§] 
and references therein. 

An alternative approach to the QCD phase 
diagram is functional methods, i.e. the func- 
tional rcnormalization group (FRG) and Dyson- 
Schwinger equations (DSEs). Naturally, these ap- 
proaches have to rely on approximations in most 
case£|, which are less controlled than those in lat- 
tice Monte Carlo simulations. These approxima- 
tions, however, are controlled by constraints such 
as symmetries and conservations laws and also by 



1 There are, however, limits where exact results are possi- 
ble, see e.g. Ref. Ql . 



comparison of the results with corresponding ones 
from lattice calculations at zero and imaginary 
chemical potential. If this procedure is followed 
carefully, one may hope to obtain meaningful re- 
sults also for the finite chemical potential region 
with Hq/T > 1. 

One of the potential advantages of the FRG and 
DSE-approach to QCD as compared with models 
like the Nambu-Jona-Lasinio model (NJL) 11], 
its Polyakov loop extended versions [12h14| and 
the Polyakov loop extended quark-meson model 
(PQM) [l5l - [l7j is the direct accessibility of the 
Yang-Mills sector. In the models, gluons are not 
active degrees of freedom and their reaction on 
the medium can thus neither be studied nor di- 
rectly taken into account. This is possible within 
the functional approaches to QCD. With func- 
tional renormalization group equations t he g luons 
of quenched QCD where studied in Ref. [13], and 
QCD at zero and imaginary chemical potential in 
Refs. [IH, H3]. With Dyson-Schwinger equations, 
lattice data for the quenched gluon where used as 
an input for the quark DSE to study chiral and 
deconfinement transitions as well as quark spectral 
functions in quenched QCD [2l|, [22| ■ These calcu- 
lations were generalized to the two-flavor case and 
to finite real chemical potential in Ref. [23| . 

In the latter study the quarks have been back- 
coupled to the Yang-Mills sector using a hard ther- 
mal loop (HTL) like approximation for the quark 
loop in the gluon DSE. In this work we will lift 
this restriction and consider the back-reaction with 
dressed quarks, as discussed in section [II] and the 
appendix. We will therefore be able for the first 
time to calculate the unquenched two plus one fla- 
vor quark and gluon propagators at finite tempera- 
ture and density. We also determine the chiral and 
deconfinement transition regions from the quark 
condensate and the dressed Polyakov loop [111 H3- 
|26| . Our results for the QCD phase diagram with 
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Nf = 2 and Nf = 2 + 1 quark flavors are pre- 
sented and discussed in section Hill In the conclud- 
ing section |IV] we summarize our results and dis- 
cuss remaining uncertainties due to our truncation 
scheme. 



II. IN-MEDIUM PROPAGATORS AND 
ORDER PARAMETERS 

The objects of interest in this work are the prop- 
agators of the quark and the gluon, from which 
we will extract their thermal properties and order 
parameters for chiral symmetry breaking and con- 
finement. In the Landau gauge and in the medium 
they are given by 

S( P ) = [i{u n + inftiCfr) + ip^A( P ) + B(p)]" 1 , 

(1) 

ZV(P) = + PlAv)^ ■ (2) 

The vector dressing functions A and C and the 
scalar dressing function B of the quark propaga- 
tor depend on the momentum p = (u> n ,p) and, 
implicitly, on temperature and chemical potential. 
The same is true for the gluon dressing functions 
Zl and Zt with longitudinal and transversal ori- 
entation with respect to the heat bath. The corre- 
sponding projectors are given by 

p%, = (i - <W) (i - M (V - Pj p L ) , (3) 

P\LV = ~ (4) 

The Matsubara modes w n are subject to antiperi- 
odic boundary conditions for fermions, lo„ = 
7rT(2n + l), and periodic ones for the gauge boson, 
Lj n = TrT2n. The quark dressing-functions will be 
determined self-consistcntly from the quark DSE, 
while the gluon dressing-functions will be given by 
a combination of lattice results and the gluon DSE 
as discussed in section Hi Al below. 

From the quark propagator, Eq. ([T]), one can 
calculate the condensate 

where Z^ is the quark wave function renormaliza- 
tion constant. For finite bare quark masses the 
resulting quantity is quadratically divergent and 
can be regularized according to 

A,, a = (^)i-— <^>«. (6) 



where the divergent part mA 2 from the strange 
quark condensate ('ipip)s cancels the correspond- 
ing part of the light quark condensate {ipip)i when 
multiplied with the ratio mi/m s of bare light to 
strange quark masses. The quark condensate is a 
strict order parameter for chiral symmetry break- 
ing in the chiral limit and serves as an indicator 
for the chiral crossover at physical quark masses. 

In [23| we presented a first calculation of the 
dressed Polyakov loop at finite real chemical po- 
tential in two-flavor QCD. This object is defined 

by mm 
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£„ = / %^ vn m^ (7) 



where (ipip)tp is the quark condensate evaluated 
at generalized, £/(l)-valucd boundary conditions 
ip(x,l/T) = e^ip(g,0) with <p £ [0,2tt[. For E±i 
to act as order parameters for deconfincment it 
is mandatory to implement the generalized, U(l)- 
valued boundary conditions only on the level of 
observables, but not in the partition function it- 
self. All closed quark loops therefore maintain the 
physical value tp = 7r, whereas tp € [0,27r[ in the 
quark-DSE. This procedure breaks the Roberge- 
Weiss symmetry, a necessary condition for the dual 
condensate to act as an order parameter for center 
symmetry breaking [2(| [27| ■ 

In the presence of a crossover for chiral symme- 
try breaking and confinement, we determine the 
maxima of the derivative with respect to the quark 
mass of the order parameters to obtain the pseudo- 
critical temperatures. 

The main difficulty in the Dyson-Schwinger 
framework is to find a truncation scheme that cor- 
rectly describes the relevant physics. In the case of 
QCD thermodynamics this means to cut off the in- 
finite tower of DSEs in such a way that the temper- 
ature and density dependence of the n-point func- 
tions that are not determined self-consistently is 
carefully approximated. Fig. Q] shows the DSE for 
the quark propagator, which depends on the fully 
dressed gluon propagator and quark-gluon vertex. 
For the quark-gluon vertex we will use a carefully 
designed expression constructed along its Slavnov- 
Taylor identity and matched to provide the correct 
renormalization group running in the quark and 
gluon DSEs. For the gluon we will use temperature 
dependent lattice data for the quenched propaga- 
tor, and unquench it by invoking the quark part of 
the gluon DSE. In particular this means that the 
gluon becomes sensitive to the chiral dynamics of 
the quark sector. 
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FIG. 1. The DSE for the quark propagator. Large 
blobs denote dressed propagators and vertices. 




FIG. 2. The full gluon DSE. Blobs denote dressed 
propagators and vertices. 
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FIG. 3. The truncated gluon DSE for N f =2 + 1 QCD. 
The yellow dot denotes the quenched propagator. (For 
interpretation of the references to color in this figure 
legend, the reader is referred to the web version of this 
Letter.) 



A. The gluon DSE 

Fig.|2]shows the untruncated gluon DSE. Studies 
of this equation at finite temperature have proven 
to be very difficult [HI, [29|, and are mostly re- 
stricted to large and infinite temperatures. An- 
other possibility to access the finite-temperature 
gluon is to use the functional renormalization 
group, where interesting first results have been 
reported recently [l8|. Although the system- 
atic errors arc still considerable, probably the 
most reliable resource, however, is lattice QCD, 
where a number of quenched studies are available 
[22l [28l . [30M32| . These errors are most pronounced 
around the deconfincmcnt transition temperature 
of the pure gauge theory, but fortunately have not 
much impact on our unquenched calculation [33| . 
Here, we will continue to use the lattice data from 
Ref. [22J as also done in our previous work [23| . 



In terms of diagrams, quenched QCD means 
that all quark loops are neglected, i.e. the last dia- 
gram of Fig.[3J One then recovers the first (second) 
order deconfincmcnt transition of SU(3) (SU(2)) 
Yang-Mills theory at the correct transition tem- 
peratures. Using the dressed Polyakov loop this 
transition has been studied in Ref. (22|, a corre- 
sponding study for the Polyakov loop potential us- 
ing FRGs is reported in Ref. [HI . In order to inves- 
tigate the phase diagram at physical quark masses, 
however, one needs to take the quark loop explic- 
itly into account. In fact, as we will see a care- 
ful treatment of the back-reaction of the quarks 
onto the gluon sector is mandatory for a mean- 
ingful study of the phase diagram. Consequently, 
this also introduces an implicit chemical-potential 
dependence of the gluon dressing functions which 
is most pronounced in their thermal masses, see 
section IIIIBI 

Technically, we implement this dependence by 
substituting the Yang-Mills part of the gluon DSE 
by the quenched propagator from the lattice, and 
merely add the quark loop. The resulting equa- 
tion is shown diagrammatically in Fig. [3l where 
we already distinguished between the (isospin- 
symmetric) light quark and the strange quark con- 
tribution. This procedure of merely adding the 
quark loops does neglect all quark loop effects in- 
side the gluon or ghost loops and is therefore only 
approximate. However, we explicitly checked that 
this approximation works well in the vacuum [35| , 
where differences to the fully self-consistent solu- 
tion are below the five percent level. This may still 
be the case at finite temperature and chemical po- 
tential. 

In [23l | we furthermore approximated the quark 
loop by taking only bare quarks into account, 
which amounts to using a hard-thermal loop 
(HTL) like approximation with only the quark- 
gluon vertex being dressed (for a review on HTL 
results see e.g. [36[). This approximation, intro- 
duced for technical reasons, is much more severe 
and can only be justified for temperatures well 
above the critical one. In this work we there- 
fore improve upon this situation and calculate the 
quark loop explicitly with the quark dressing func- 
tions obtained from the quark-DSE. By doing so, 
the quark and gluon DSE become coupled and have 
to be solved simultaneously. This means also, that 
the gluon becomes sensitive to the chiral dynamics 
in the quark sector and in particular to the chiral 
transition. The coupled system of DSEs in Figs. [T] 
and Fig. [3] can be written as 
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[s'GO] 1 = z{ [s[( P )\ 1 + c F ziz{ F j£g^sf(i) g rl(P,p 2 ,q 2 )D^(q) 7 
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where g = p — I, is the quark propagator for 
one flavor / = u, <i, s, Cf = f is the Casimir oper- 
ator, r„ the dressed quark-gluon vertex, Z\p and 
Z2 are the vertex and wave function renormaliza- 
tion constants, and the Matsubara sum as well as 
the integration over the loop three-momentum I is 
represented by £ ; = TJ2 n J j§pp- For the mo- 
menta we often abbreviate gencrically p = (p, ui p ) 
and p 2 = p 2 + to 2 . With the transversal and longi- 
tudinal projectors from Eqs. (|3l4p . the quark loop 
can be decomposed as 



(8) 



The medium effects in the gluon-DSE manifest 
themselves predominantly in contributions to the 
thermal mass of the gluon. As a result, the quark 
loop can be split into a finite part that is similar to 
the vacuum quark loop and an infrared divergent 
thermal part that is proportional to 1/p 2 : 



n T > L (f 



( m th L ) 



'-reg 



(9) 



I 

with the infra-red regular part H reg . and thermal 
masses 



(10) 



A technical problem associated with the use of a 
sharp momentum cut-off A in the quark loop is the 
appearance of quadratic divergencies. This prob- 
lem is particularly troublesome in the medium, 
since the A 2 term appears like a squared ther- 
mal mass and has to be carefully removed with- 
out spoiling the T 2 and /i 2 terms which carry the 
physics. We do this using the Brown-Pennington 
projection method [37| . as detailed in Appendix |A"1 
After removing the quadratic divergences and pro- 
jecting on longitudinal and transversal parts the 
quark loop in the medium is given by 



1 



D q (l)D q (q) 



A{l)A{q) 



r s (l,q) U-^-lq) +T 4 (l,q)l-q 



B(l)B(q) [r 4 (/, q) - T s (l, q)] + C{l)C{q) [-W 2 (T 3 (l, q) + T 4 (l, <?))] 



(11) 



(12) 



with q =p-l and D q (p) = p 2 A 2 {p) + w 2 C 2 (p) + 
B 2 (p) for one quark flavor and with the abbre- 
viation Co n = uj n + ifJ-- The quark-gluon vertex 
has been split in three scalar functions, T, T s and 
T 4 which will be explained below. For the limit 
p — >• we find the thermal masses. Of course, the 
expressions (fTU)) together with (fTTj) , (fT"2"|) also have 
the correct high-temperature limit, where the HTL 



results (mf h f = and «.) 2 = g (t 2 + 3fT), 
for one flavor, are recovered. 



Some technical details of the evaluation of this 
object are given in App. [X] It is noteworthy that 
the quark loop docs not contribute to the transver- 
sal thermal masses, in contrast to the Yang-Mills 
part of the gluon self-energy. 
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B. The quark-gluon vertex 

An important quantity that appears in both, 
the quark and gluon DSEs, is the dressed quark- 
gluon vertex. At zero temperature the nonpertur- 
bativc structure of this vertex has been explored 
to some degree in the past (38l - l40j . however not 
much is known about the behavior at finite tem- 
perature and density. Therefore in order to close 
the system of DSEs we have to resort to an edu- 
cated guess for the form of the quark-gluon inter- 
action. Fortunately, there are constraints which 
serve as important guides in the construction of 
such a guess. On the one hand, at large momenta 
temperature and density effects are exponentially 
suppressed and one can entirely rely on the vacuum 
structure of the vertex. It is then well known, that 
the vertex has to combine with the gluon dressing 



functions such that the resulting combination of 
dressing functions runs like the running coupling in 
the ultraviolet. In addition, the vertex has to sat- 
isfy its Slavnov- Taylor identity. This identity has 
not been solved yet in QCD, however important 
guidance can be obtained from the correspond- 
ing Abclian Ward identity, which serves to express 
part of the vertex dressing in terms of the quark 
dressing functions. All these constraints have al- 
ready been taken into account in previous works 
[2l| - [23| . We therefore employ the same construc- 
tion which uses the first term of the Ball-Chiu ver- 
tex, satisfying the Abelian Ward-Takahashi iden- 
tity, multiplied with an infrared enhanced func- 
tion r(p 2 , k 2 ,q 2 ) that accounts for the non- Abelian 
dressing effects and the correct ultraviolet running 
of the vertex. The resulting expression reads 



IV(p, k;q) = 7 M ■T(p 2 ,k 2 ,q 2 ) ■ (8^4 



C(p) + C(q) A(p)+A(q) 
o + <W o 



rV,fcW) = 



di x f (3 a{pi)\n[x/A 2 + 1] 

d 2 + x A 2 + x \ Att 
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(13) 



where p and k are the fermion momenta, q is the 
gluon momentum and d\ as well as the scales d 2 
and A are parameters. Whereas d 2 and A control 
the renormalization group running of the vertex 
function from the large into the low momentum 
region, d\ controls the strength of the quark-gluon 
interaction at small momenta and therefore the 
amount of quark mass generation in the hadronic 
phase. In the ultra-violet 8 = — ^ ^ iN N ^ SN is 
the anomalous dimension of the vertex and ft® — 
11N "~ 2Nf . In Eqs. (|11I12[) we have used the ab- 
breviations T s (p,q) = ^M±fKgJ and T^(p,q) = 

C(p)+C(q) 

1 

The squared momentum variable x is chosen to 
be q 2 in the quark DSE, but p 2 + k 2 in the quark 
loop. This change in the momentum dependence 
is well justified by the need to maintain multiplica- 
tive renormalizability in the gluon-DSE [35[ . Note 
that, since the Ball-Chiu part is determined from 
the quark dressing functions our vertex does in- 
clude effects from temperature and chemical po- 
tential. 



C. Strange quarks and unquenching effects 

The evaluation of the quark loop with dressed 
quarks allows us to add strange quarks in a non- 
trivial way. In general the light and strange quarks 



couple via unquenching effects in QCD, as shown 
in Fig. [3l In the vacuum, results for the propa- 
gators of Nr = 2 + 1 QCD have been discussed 
in Ref. [35|. Here we present the first finite tem- 
perature and chemical potential study of the prop- 
agators in the unquenched theory. However, we 
also wish to emphasize that not all the unquench- 
ing effects arc included yet. In addition to the 
ones in the gluon propagator there are also cor- 
responding effects in the quark-gluon vertex. Di- 
agrammatically, part of these can be written as 
hadronic contributions like pion and kaon exchange 
[39^ In Ref. [4l[ it has been argued that these ef- 
fects are necessary to obtain critical scaling beyond 
mean field at the chiral phase transition temper- 
ature in the chiral limit. In our case of physical 
quark masses, however, we expect their influence 
on the critical temperatures to be small. Further- 
more, since mesons have zero quark number the 
chemical potential dependence of the meson ex- 
change is expected to be small as well. In this 
work we therefore do not include these contribu- 
tions explicitly, but rather absorb their effects in 
our vertex ansatz (fT3")) . While this strategy may 
work rather well for meson contributions, it pre- 
sumably fails to account for diquark and baryon 
contributions at large densities and low temper- 
atures. This is the region of the phase diagram 
where nuclear effects are expected to play a ma- 
jor role [H and results in our present truncation 
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scheme should be interpreted with great care. We 
will come back to this point below. 

In our calculation we only vary the chemical po- 
tential of the up/down quarks and set the strange 
quark chemical potential to [i s = 0, since the net 
number of strange quarks in heavy ion collisions 
vanishes. The light quark mass can be fixed by cal- 
culating the pion mass via the Gcll-Mann-Oakes- 
Renncr relation and the Pagel-Stokcr approxima- 
tion for the pion decay constant. This gives a mass 
of m u / c i k, 2 MeV at a renormalization point of 80 
GeV. We furthermore fix the strange quark mass 
by its ratio to the light quark mass by choosing 
™° = 27. For the vertex parameters we take 

di = 7.5 GeV 2 , d 2 = 0.5 Go V 2 and A = 1.4 GeV. 
For the coupling constant we take a(fi) = 0.3. 
With this set of parameters and with two mass- 
less and one strange quark, the pion decay con- 
stant in the vacuum is f v m 88 MeV and at the 
physical point of quark masses at vanishing den- 
sity, T c k 156 MeV. Those numbers agree very 
well with expectations, which is in contrast to sim- 
pler truncation schemes where usually the critical 
temperature comes out rather small [42|. 

Note again, however, that unquenching effects 
in the vertex are taken into account indirectly via 
the quark dressing functions A and C that appear 
explicitly in the vertex construction (fT3"]) . We em- 
phasize that it is nontrivial that this vertex con- 
struction leads to correct quenched [22j as well as 
unquenched transition temperatures as will be de- 
tailed in the next section. 



III. RESULTS 

Condensates, thermal gluon mass and 
phase diagram for Nf — 2 
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FIG. 4. Quark condensate, dressed Polyakov loop 
(dual condensate) and thermal electric gluon masses 
for vanishing quark chemical potential with Nf — 2. 
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FIG. 5. The phase diagram for Nf = 2 in the temper- 
ature and (up/down) quark chemical potential plane. 
Shown are the transition lines determined using the im- 
proved quark loop of this work (see legend) compared 
with our previous results (light colors in background) 
in the HTL-approximation of the quark loop [23f [ . 



We first discuss the effect of coupling the two- 
flavor quark and gluon DSEs together as explained 
above. Here we will encounter some qualitative 
effects, which arc still valid in the Nf = 2 + 1 case 
discussed below. 

The quark condensate and dressed Polyakov 
loop which are shown in Fig. [J] for two light quark 
flavors at vanishing chemical potential exhibit a 
crossover behavior around T c « 200 MeV, with the 
transition region for chiral symmetry restoration 
and deconfinement being almost equal. We clearly 
identify the two most important effects from in- 
cluding the quark loop: a reduction of T c and 
a change from a first order phase transition to a 
crossover. Fig. 2] also shows the contribution from 
the quark loop to the electric part of the thermal 
gluon-mass as defined in Eqs. (fTUj) . (fT2|) . The ther- 
mal mass is small in the hadronic phase, shows 
a sharp rise in a temperature range around the 



transition temperature and converges to the T 2 - 
bchavior expected from HTL-calculations in the 
high temperature phase. This behavior can be ex- 
pected, since the quark loop has an inverse depen- 
dence on the quark mass-function. It is therefore 
suppressed by dynamical chiral symmetry break- 
ing in the hadronic phase and becomes large in 
the high temperature phase, where the quark con- 
densate becomes small. As compared to a trunca- 
tion scheme where the quarks in the loop are not 
dressed, i.e. the HTL-likc approximation we used 
in [23| , the quark condensate shows a steeper tran- 
sition around T c . This behavior originates from the 
growing thermal gluon mass due to the quark loop, 
which in turn decreases the interaction strength in 
the quark DSE. Therefore, the back coupling of the 
quark and gluon DSEs accelerates the chiral phase 
transition. 

In Fig. [5] we present our results for the chiral 
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and deconfincmcnt transition lines determined in 
this work and compared to our previous results of 
Ref. [23| . where we used the HTL approximation 
for the quark loop. Note that the deconfinement 
transition is determined solely from since 
the corresponding calculation of £+1 is plagued 
by considerable numerical uncertainties we refrain 
from showing the results. We first observe that 
the pseudo-critical temperatures for restoration of 
chiral symmetry and deconfinement are now closer 
together, in comparison to the HTL results where 
a difference of up to 10 MeV was observed. This 
closer matching of the critical temperatures can be 
explained by the gluon becoming sensitive to the 
chiral phase transition due to the explicit appear- 
ance of the quark mass functions in the full quark 
loop calculation. Moreover, the critical end-point 
has moved to a smaller quark chemical potential 
and a larger temperature. We attribute this to the 
accelerated chiral transition due to the improved 
back-coupling of the quarks onto the gluon, as ex- 
plained above. The effect is substantial and under- 
lines once more the general importance of details in 
the back-coupling effects from the quarks onto the 
Yang-Mills sector of QCD. As a result we obtain a 
ratio of critical quark chemical potential to critical 
temperature at the end point of fj, q EP /TEP ~ 1.1. 
Thus for Nf = 2 QCD this ratio stays above one 
(although not much) in agreement with the claim 
made in Ref. (23|. We will see below that the ad- 
dition of the strange quark further increases this 
ratio. 



B. Condensates and phase diagram of 

Nf = 2 + 1 QCD 

In addition to the two light quark flavors we now 
also add the heavier strange quark to the system. 
Recall, however, that following the conditions in 
a heavy ion collision we do not include a chemi- 
cal potential for the strange quark in our calcula- 
tions. We discuss the effect of this choice below. 
Let us first have a look at the resulting unquenched 
gluon propagator. The main effect of the tempera- 
ture and chemical potential dependent quark loop 
onto the gluon is the change of the thermal elec- 
tric mass. We plotted this quantity, normalized by 
its asymptotic value, as a function of temperature 
and up/down-quark chemical potential in Fig. |BJ 
One clearly sees the continuous change for small 
chemical potential similar to the two-flavor case in 
Fig. |U For larger chemical potential the transi- 
tion becomes steeper until it becomes discontinu- 
ous in the vicinity of the critical endpoint of the 
chiral transition, discussed below. This behavior 
underlines again the fact that the chiral transi- 
tion in the quark sector of the unquenched theory 
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FIG. 6. The thermal gluon mass normalized by its 
asymptotic (HTL) value (omitting g 2 ) as a function of 
temperature and chemical potential. 




T [MeV] 

FIG. 7. The regularized condensate compared to lat- 
tice results from [l[. 



also drives the changes in the Yang-Mills sector of 
QCD. In the asymptotic region at large temper- 
ature and/or chemical potential we finally recover 
the HTL-result m 2 th ~ T 2 + 3fj 2 /n 2 for the thermal 
gluon mass. 

Let us now focus on the quark sector of QCD and 
come back to the behavior of the quark condensate 
at zero chemical potential. Since we have two plus 
one quark flavors at our disposal we can directly 
compare with corresponding lattice data at van- 
ishing chemical potential. In Fig. we compare 
the renormalized quark condensate A; jS defined in 
Eq. ([6]) with lattice data from Ref. [l[ . The agree- 
ment is good except for the region above the phase 
transition, where the DSE results are larger than 
the lattice data. We attribute this discrepancy to 
the truncation of our quark-gluon vertex, which 
employs medium effects only in the Ball-Chiu part, 
while the ansatz function T is temperature inde- 
pendent. This is in contrast to [43j , where it has 
been suggested that the parameter d\ of the vertex 
infra-red part is reduced above T c . Therefore, ne- 
glecting the T-dependence of the vertex infra-red 
part leads to an over-estimation of the interaction 
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FIG. 8. The light (lower surface) and strange (upper 
surface) quark condensate as a function of temperature 
and chemical potential. 



strength in the chirally symmetric phase, and thus 
to too large condensates. Also the scalar parts in 
the vertex are neglected, which are known to be 
important when it comes to the details of dynami- 
cal chiral symmetry breaking and restoration [38j . 
Additionally, the (neglected) temperature depen- 
dence of ex plic it hadronic parts of the vertex may 
play a role [4l[. Nevertheless, we find our present 
truncation sufficiently accurate in particular with 
respect to the steepness of the crossover to justify 
its extension into the finite chemical potential do- 
main. 

In Fig. [5] we show our results for the light and 
strange quark condensates in the temperature and 
chemical-potential plane. The condensates are not 
renormalized and shown on different scales for the 
sake of clarity. In the background of the figure 
we see again the crossover behavior of the light 
quark condensate at small values of the chemical 
potential. This turns into a first order transition at 
large chemical potential, visible in the foreground 
of the figure. The strange quark condensate shows 
a behavior similar to the light quark condensate 
around the phase transition, but continues to melt 
at larger temperatures. Certainly, this behavior 
is caused by the larger explicit breaking of chi- 
ral symmetry for the strange quark. At densities 
where the light quark shows a first order phase 
transition, the strange quark also shows a discon- 
tinuity at the same temperature. This is inherited 
from the discontinuity in the gluon propagator, see 
Fig. |51 and demonstrates how the gluon couples 
light and strange quarks. 

In Fig. [9] we show the resulting phase diagram 
for 2 + 1 flavor QCD in our approximation scheme. 
As compared to the two-flavor case, Fig. [SJ we find 
a chiral critical endpoint at slightly larger values of 
chemical potential. Thus, although the additional 
strange quark lowers the transition temperatures 
by roughly 45 MeV, the melting of the condensate 
is slightly less steep as for the Nf = 2 case. As 
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Chiral 1 st order 
CEP 
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FIG. 9. The phase diagram for two plus one flavors. 
The light colors (top lines) show the Nf = 2 results as 
a comparison. 



a consequence, one hits the first order transition 
region at larger values of the chemical potential. 

At (i = we furthermore observe a slightly wider 
spread between the chiral and the deconfinement 
transition for the Nf = 2 + 1 case as compared to 
Nf = 2. This effect may be explained by the dif- 
ferent strength of the back-reaction of the quarks 
onto the gluon sector. For Nf = 2+1 the back- 
reaction is stronger, resulting in a wider crossover 
transition region which in turn allows for a wider 
spread between the transition temperatures. On 
the other hand, the quenched calculation with no 
back-reaction is entirely driven by deconfinement 
and both transitions happen at the same temper- 
atures [22| . 

So far, we set the strange quark chemical po- 
tential to zero. Setting it equal to the light quark 
chemical potential has very little effect on the re- 
sults. 

From the chiral transition line one can extract 
the curvature k, which is defined by 



T C ( M ) = T c (0) 



1 - K 



Tc(0) 



(14) 



as the coefficient of a power expansion around /.t = 
0. This quantity can be measured on the lattice 
despite the fermion sign problem, and therefore 
serves as a comparison of the chemical potential 
dependence of the transition line. Two different 
lattice groups determined the curvature recently 
and found k « 0.059(2) (4) | and k k 0.059(20) 
Q for 2 + 1 flavors. In comparison with the Nf = 2 
result k w 0.051(3) of Ref. Q, there is a general 
trend of increasing curvature when more flavors are 
taken into account in agreement with a large N c - 
analysis [45[. In our case, we extract k by fitting 
Eq. (jT4")) to the chiral transition line in the region 
\i e [0, 25] MeV and n E [0, 50] MeV for the case 
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N f 


CEP 


t c (m = 0) 


K 


2 (HTL) 
2 

2+1 


(280,90) 
(171,154) 
(190,100) 


183 MeV 
202 MeV 
156 MeV 


0.23 
0.41 
0.30 



TABLE I. Location of CEP and the curvature for Nf = 
2 and Nf = 2 + 1 flavors as well as the Nf = 2 flavor 
result in the HTL approximation of Ref. |23l j. 

of the HTL-like quark loop. We also checked that 
the order (/la/T) 4 is negligible in these intervals. 

Our results for the curvature k from the HTL- 
truncation of Ref. [23[ and the two results of this 
work are summarized in Tab. U together with our 
findings for the location of the critical end-point 
(CEP). We find that the back-coupling increases 
the curvature compared to the HTL-like trunca- 
tion. This is easily understood from the shape of 
the transition lines of Fig. [51 where the CEP from 
the fully back-coupled results is almost on top of 
the transition line of the HTL-approximated re- 
sult, while at the same time the transition tem- 
perature at /z = is significantly larger. In gen- 
eral, our values for k agree in size with QM and 
PQM calculations [46|,|47| but are much larger than 
the lattice results. Part of this discrepancy may 
be attributed to volume effects in the lattice cal- 
culations HU, another part of it may be due to 
limitations of our present approximation scheme. 
This is particularly apparent in the A'j-dependence 
of the curvature. As discussed above, the lattice 
results indicate a significant increase of curvature 
with Nf, whereas our results indicate the opposite. 
The reason for this discrepancy may be found in 
the incomplete treatment of unquenching effects in 
our present approach. While we took into account 
quark loop contributions in the gluon polarization, 
the inclusion of corresponding effects in the quark- 
gluon vertex is numerically demanding and will be 
investigated in future work. 



IV. CONCLUSION AND OUTLOOK 

In this work we have determined the phase 
boundaries of Nf = 2 + 1 QCD in the temperature 
and (up/down) quark chemical potential plane. 
We have solved a coupled set of Dyson-Schwinger 
equations for the Landau gauge quark and gluon 
propagators. For the gluon we worked with lattice 
results for the temperature dependence of the elec- 
tric and magnetic part of the quenched propaga- 
tor and added quark loop effects determined in the 
DSE-framework. For the quark-gluon interaction 
we worked with a construction which takes into ac- 
count the correct renormalization group running of 
the vertex combined with leading order tempera- 
ture and chemical potential effects as encoded in 



its Slavnov- Taylor identity. 

From the temperature and chemical potential 
dependence of the quark propagator we extracted 
information on the chiral as well as the deconfinc- 
ment transition encoded in the quark condensate 
and the dressed Polyakov loop. As a result we 
find a crossover behavior at small chemical poten- 
tial in quantitative agreement with corresponding 
lattice results. We also find a critical endpoint 
at (T EP ,n EP ) « (130,145) MeV followed by a 
first order transition. We therefore confirm our 
previous claim of the nonexistence of a CEP for 
fx EP /T EP < 1 |23j in agreement with extrapolated 
results of lattice QCD [J, [f|. 

The calculations presented in this work are, to 
our knowledge, the first effort based on continuum 
QCD to determine the phase diagram with realistic 
up/down and strange quark masses combined with 
a dynamical gluon sector. In this respect we wish 
to emphasize again the importance of the back- 
coupling effects of the quarks onto the Yang-Mills 
sector of the theory. These effects are responsible 
for the change of the order of the transition from 
first order in the quenched theory to the crossover 
behavior of full QCD and the corresponding reduc- 
tion of the (pseudo-) critical temperatures. Both 
phenomena are nicely reproduced in our present 
truncation scheme. 

Nevertheless, this is not the end of the story and 
from the present approximation scheme no firm 
conclusions on the precise location of the CEP nor 
on its very existence can be drawn. The reason is 
to be found in our present expression for the quark- 
gluon interaction. Although it does take into ac- 
count the above-mentioned constraints from QCD 
itself, it lacks in detail when it comes to the precise 
dependence on temperature and chemical poten- 
tial encoded in hadron back-reaction effects onto 
the quark propagator, as discussed in section iHBl 
In the DSE-framework, these have only just begun 
to be explored [4l|. In particular, baryon effects 
may be important when it comes to large chemical 
potential, as has been emphasized in Q and explic- 
itly checked in a two-color framework in Ref. [48j. 
Furthermore, at large chemical potential inhomo- 
geneous phases may be energetically favored over 
homogeneous ones, see e.g. (49l . [50j and references 
therein. Whether these effects show up in addi- 
tion to the existence of a CEP or instead of the 
CEP is an important open question that needs to 
be investigated further in the future. 
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Appendix A: Details of the quark loop 
calculation 

1. Vacuum 

The renormalization of the quark loop that is 
necessary to remove the logarithmic divergence is 
done in the vacuum, since medium effects do not 
lead to new divergences. We require that at the 
renormalization point £, the self energy vanishes: 



n ren (p 2 ) = n(p 2 )-n(c 2 ), (Ai) 

where we choose £ = 10 GeV. 

When solving the system of coupled DSEs in 
the vacuum, one notes that with the normal kind 
of vertex ansatz which only depends on the gluon 
momentum, as it is used in rainbow-ladder trunca- 
tion, a function which allows for dynamical chiral 
symmetry breaking is hard to find. The reason for 
this is that the quark loop becomes very strong, 
as it is multiplied with the vertex dressing func- 
tion which only depends on the external momen- 
tum here. Using a vertex which depends on all 
momenta, which a realistic vertex certainly would, 
results in a weaker quark loop and therefore allows 
for chiral symmetry breaking. Also, multiplicative 



renormalizability is spoiled by a dependence of the 
vertex on the external momentum. 



2. Medium 

We already mentioned above that the dressing 
functions in the gluon DSE are plagued by spu- 
rious quadratic divergences which appear when a 
hard cut-off is used. In the vacuum they can be 
subtracted in several ways, in the medium one has 
to be careful since they appear like thermal masses: 

n L (0) = a-A 2 + &-T 2 + c-/i 2 , (A2) 

and one has to be careful to only remove the first 
term. We solve this problem by using a Brown- 
Pcnnington projection 

rWp) (re9 - } = n^(p) - v^n a/3 ( P ), (A3) 

which has been introduced in [37|. With the 
Brown-Pennington projector it is vital to use an 
0(4)-invariant cut-off, or the quadratic divergence 
appears again. We achieve this by performing the 
Matsubara sum and the integral as follows: 



n n =-N-l *{ ^ 

(A4) 

where A' = ttT(2N + 2) is a temperature- 
dependent cut-off, that is chosen close to the fixed 
cut-off A by taking N from the Matsubara mode 
closest to A. 
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